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ABSTRACT 

The observation of energetic X-ray emission from black holes, inconsistent with thermal emission from an 
accretion disk, has long indicated the presence of a "corona" around these objects. However, our knowledge 
of the geometry, composition, and processes within black hole coronae is severely lacking. Basic questions 
regarding their size and location are still a topic of debate. In this letter, we show that for black holes with 
luminosities L > 10~^ LEdd - characteristic of many Seyferts, quasars, and stellar-mass black holes (in their 
brighter states) - advanced imaging and timing data strongly favor X-ray emitting regions that are highly 
compact, and only a few Gravitational radii above the accretion disk. The inclusion of a large number of 
possible systematics uncertainties does not significantly change this conclusion with our results still suggesting 
emission from within ^ 20rg in all cases. This result favors coronal models wherein most of the hard X-ray 
emission derives from magnetic reconnection in the innermost disk and/or from processes in the compact base 
of a central, relativistic jet. 



1. INTRODUCTION 

Observations of accreting black holes (BHs) often indicate 
the presence of hard X-rays having energies much greater than 
the expected thermal peak of the accretion disk. This is true 



not on ly for Active Galactic Nuclei (AGNs; e.g. Elvi s et al. 
T978]l but also for X-ray binaries (XRBs; e.g. White & Holt 
IWl^. 



It has been widely postulated that hard X-rays are the prod- 
uct of inverse Compton scattering of seed photons from accre - 
tion disks by hot "coronae" (e.g. [Haard t & Maraschi 1991), 
however despite our observational (e.g. 'Zdzia rski et al._1999_ ) 
and theoretical (e.g. Schnittman et al. 2012, and references 
therein) efforts, there are still a number of fundamental ques- 
tions regarding the geometry, or composition of accretion- 
disk coronae. 

At low luminosities {L < lO^'^Lsdd^ where L^dd is the 
Eddington luminosity), the hard-X-rays are thought to origi- 
nate from a hot, quasi-spherical and large (^ IQOGM/c^ ^ 
lOOrg) advection-dominated accretion flow (A DAF) which 
replaces the innermost disk around the BH (e.g. |Narayan"&| 



Yi|1994| l. On the other hand, at luminosities L > lO'^Lsdd, 
observations of XRBs and Seyfert AGNs show the clear pres- 
enc e of a more comp act X-ray emitting source (the corona; 
e.g. [McHardy et al.|20 06) on top of a geometrically -thin, op- 
tically thick accretion disk extending as far in as the radius of 
the innermost stable circular orbit (ISCO). 

Amongst the clearest evidence for the co-existence of the 
corona with the disk is the presence of Fe K-shell emission 
lines together with oth er reflection f eatures in the X-ray spec- 
tra of these sources ( |Miller|[2007[ and references therein). 
From the profile of these emission lines we are able to obtain 
information on the effects of strong Doppler and gravitational 
redshifts natura l to the spa ce-time around BHs and measure 
their spin (see Miller|2007 for a review). 

Having a tool to probe the innermost regions around BHs 
and to measure their spin (a — cJ /GAP, — 1 < a < 1) is 
of particular importance as the distribution of spin in nearby 
supermassive black holes (SMBHs) can be used to distinguish 
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betwe en various BH growth scenarios (e.g.|Berti & Volonteri 
2008|l, and in the case of stellar mass BHs, their spin distri- 



bution relates to t he distribution of an gular momentum of the 
progenitors (e.g. Miller et al.||201l| l. Knowledge of spin is 
also essential if one wishes to test if relativistic jets can be 
powered by the spin energy of BHs via th e Blandford-Znajek 
process (e.g. Narayan & McClintock 2012) or, in the near fu- 
ture, if t here is evidence for violatio ns of the No-Hair Theo- 
rem (e.g. Johannsen & Psaltis|201 1[ ). 

The location and compactness of the corona - which is irra- 
diating the disk and giving rise to the reflection features - are 
thus of paramount importance in the study of strong gravity. 
Nonetheless, in theoretical consideration of the corona it is 
often assumed to be compa ct and located in the inner regions, 
< lOrg around the BH (e.g. |Miniutti & Fabian|2004[ and ref- 
erences therein). This assumption has consistently been suc- 
cessful in describing the spectral and timing properties of BHs 
at the highest signal-to-noise available, however strong devi- 
ation from th is could affect our ability to probe strong gravity 
around BHs Pauser et al.|2013) . 

Furthermore, coronae are know n to play an imp ortant role 
in state transitions in XRBs (e.g. 'Reis et al."2013i and ref- 
erences therein); and jets in both XRBs and AGN s have of- 
ten be en strongly linked with the corona (e.g. Merlon i et al. 
|2003|l. If the corona is the base of such a jet (Falcke & Bier- 
mann 1999), then the amount and the form in which energy 
escapes through the jet is largely dependent on the enigmatic 
properties of the corona, highlighting the importance in estab- 
lishing the size and location of these regions. 

In this letter, we explore the results aimed at measuring 
the physical sizes and locations of coronae around luminous 
accreting SMBHs (L > lO^^LEdd typical of bright nearby 
Seyferts and quasars). We find that in all cases the coronae 
are located significantly closer than ^ lOrg from the accre- 
tion disks and their sizes are suggestive of being highly com- 
pact. These results give credence to the assumptions often em- 
ployed in most reflection-based studies of BHBs and AGNs. 

2. LOCATION BASED ON REVERBERATION LAGS 

With the discovery of a soft X-ray reverberation lag of 
- 30s in the AGN IH 0707-495 ( |Fabian et al.|2009] l, a new 
method to measure the distance (-Dcor) between the X-ray 
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Fig. 1. — Left: Representative probability distribution for the distance 
between the accretion disk and the corona (Dcor) based on soft X-ray rever- 
beration lags. Right: Probability distribution for the combined sample of 17 
Seyferts (blue). The tail extending to ~ 15rg comprises of RE J1034-I-396, 
NGC 4151 and ESO 113-G()10. Shown in red is the distribution for sources 
with lag detection at > 3cr, with the exception of NGC 4151. 

emitting region and the accretion disks was established. The 
decomposition of the spectrum of IH 0707—495 into direct 
and reflected emission components supported interpreting the 
observed lag as a signature of the light travel time between 
the corona (the hard X-ray continuum source) and the accre- 
tion disk (which reprocess the flux into a reflection component 
with atomic spectral features). 

Despite the success that this "reflection" model had in ex- 
plaining the timing and spectral features of IH 0707-495 
( Fabian et al.||2009[ l, a further interpretation was soon made 
available where the lag was due to the reverberation caused 
by scattering of X-rays passing through an absorbing medium 
that happen to lie in the line of sight to the central AGN 
( [Miller et al.||20l6) l. However, the subsequent discovery of 
similar lags in a large number of sources (see below), span- 
ning a large range of masses, Eddington fractions, and incli- 
nation angles, is likely incompatible with the specialized par- 
tial covering geometry required to explain the lags. Thus, the 
original light-crossing argument remains the best interpreta- 
tion for these soft lags. For a detailed description of soft X 
ray reverberation techniques and its a pplications see Zoghbi 
|eraL] ( |20T2l l; |Wilkins & Fabian] ( |20T3] l. 

At the time of writing, there were a total of 17 radio-quiet 
Seyferta^with measured soft lags that met our luminosity cri- 
teria. This sample is shown in Table 1 together with estimates 
for the mass and luminosity as a fraction of Eddington to each 
source collected from the literatur e. The lag s compiled here 
are taken from the work o f De M arco et al., (2012); Zoghbi 
eraL](|20T2l i ; [Fabian et al| ( |2012p ; irackett et al.| ( |2012j i and 



Kara et al.[ (|2013|l, whereas the masses a re primarily from 
^ou et al. (2010); Zhou & Wang (2005); B entz et al.[d2006 i 
and O nken et al., ( |2007) . As Seyferts are known to be highly 
variable, the luminosities shown here are meant to be seen 

^ The LLAGN NGC 4395 is not used in this work as it accretes constantly 
at L/L^dd ~ 1.2 X 10""^ and is possibly already in the ADAF regime and 
thus intrinsically different from the more luminous systems. 



Fig. 2. — BH mass (top) and L/L^jdd (bottom) as a function of disk-corona 
distance. The dashed line marks the ISCO of a non-rotating BH (6rg). 14/17 
sources in our sample suggest a corona close to the accretion disk. These 
results are independent of the mass or m as long as the system are accreting 
above a few percent of Eddington. The three singled out sources corresponds 
to those peaking at ~ 9rg (§ 2) and the sources with lag detection at > So- 
level are shown with square markers. 

only as a rough indication of the average luminosity of each 
system. 

Due to the scale-invariant nature of the physics governing 
accretion flows onto BHs, the lag time scale t depend primar- 
ily on the mass of the central BHj t ^ GM/c' = Vg/c. 
We estimate the probability distribution for the distances be- 
tween the coronae and accretion disks by drawing samples 
from the probability distribution of the lag times and masses 
as obtained from the literature and reported in Table 1 . For 
NGC 4051, MRK 766, RE'J1034h-396, NGC 7469, Mrk 335, 
PG 1211H-143, NGC 3516, NGC 5548 and NGC 4151, their 
masses were assumed to have a Gaussian distribution with 
standard deviation given by the average of the positive and 
negative errors shown in Table 1, when the reported errors 
were asymmetric. For the remaining sources with the excep- 
tion of ESO-1 13-GOlO, the logarithm of the masses were as- 
sumed to have a Gaussian distribution with the symmetric er- 
rors reported in each reference. Finally, for ESO-1 13-GOlO, 
we assumed a uniform distribution in the range given in Ta- 
ble 1. 

An example of a derived probability distribution based on 
10^ samples is shown in Fig. 1 (left) for NGC 405 1 . We report 
on the 5th column of Table 1 the median of Dcor with the la 
error taken to be the standard deviation of each source sam- 
ple. Figure 1 (right) shows the distribution for the combined 
sample, where it becomes clear that there is a peaks at ^ 2rg 
and a further, minor peak at ^ 9rg with a tail that extends up 
to ^ 15rg. 

The 3 sources giving rise to the second peak at larger 
Dcor are NGC 4151, RE J1034H-396 and ESO 113-GOlO. 
There is currently uncertainty on the mass of the SMBH in 

^ This simple formula does not take into account the effect of the Shapiro 
time delay; however, this correction is of the order of 0.5rg. The appar- 
ent distance between the corona and reflector, calculated using the simplified 
formula and neglecting the Shapiro delay, is therefore larger than the actual 
distance I Wilkins & Fabian 2013 for details). 



Table 1: Summary of sources used in this work. Tiie table is divided into two with the top half referring to the 17 Seyferts and the bottom half 
to the 5 microlensed quasars. The various columns are (1) Object name; (2) black hole masses in Solar units. For the quasars, their masses are 
in logarithmic scale as quoted in the literature; (3) representative Eddington fraction; (4) lag time in seconds (top) or logarithmic of the size of 
the X-ray emitting region in centimetres (bottom); (5) The median of the distance, < Dcor >, between the corona and the accretion disk in 
gravitational units (top) or the median of the size of the corona < 5* > also in gravitational units (bottom). 



Object 
(1) 



Mass (X 10'' Mo) 
(2) 



L/LEdd 

(3) 



Soft Lag Time (s) 

(4) 



NGC4051 


-, 70-1-0.55 

^•'■J-0.52 


Mrk 766 


1 7fi+l-56 


Ark 564 


1.86+tfr 


MCG-6-30-15 


2.olS:| 


IH 0707-495 


9 q-l-5.1 


RE 11034-^396 


4tl 


NGC 7469 


12.2 ±1.4 


Mrk 335 


26 ±8 


PG 1211-1-143 


40.5t?2'i 


NGC 3516 


31.7+4'^ 


Mrk 841 


75.91111 


NGC 5548 


44.2tlis 


NCC 6860 


39+^? 


Mrk 1040 


43.7l?-.^5 


NGC 4151 


45.71^7 


IRAS 13224-3809 


5.81^.^5 


ESO 113-GOlO 


4-10 



0.03 
0.55 

1 
0.15 

1 

1 
0.32 
0.36 
0.28 
0.013 
0.44 
0.034 
0.08 
0.04 
0.04 

1 

1 



16 ±8 

43 ±18 

26 ±15 

22 ±11 

24 ±11 

150 ± 90 

130 ± 60 

210 ±100 

250 ± 90 

380 ± 170 

360 ± 130 

260 ±115 

145 ± 72 

380 ± 170 

2100 ± 500 

92 ±31 

325 ± 89 



Object 
(1) 



\og(MBH/MQ) 

(2) 



L/Lsdd 

(3) 



Log(Size/cm) 
(4) 



Q 2237-^0305 
RXJ1131-1231 
QJ0158-4325 
HE 1104-1805 
HE 0435-1223 
PG 1115-1-080 



8.68 ±0.36 
8.32 ±0.62 

8.2 ±0.2 
9.37 ±0.33 
8.76 ±0.44 

9.1 ±0.2 



0.44 
0.03 
0.4 
0.36 
0.11 
0.37 



+0.34 

-0.29 



15.46 
14.04-14.68 

14 3+"-'' 

J^*-J-0.5 

14.2-15.0 
< 15.07 



< -Dcor > (rg) 

(5) 



1+0.6 
-0.5 

3.8t?:«3 

2Atr2 
9 1-1-0.9 
^■J^-0.7 
-1.5 
-0.9 
5-1-3.4 

— 2 5 

2.2 ±6.5 
1 p+0.5 



1.9;! 



1.9^ 
6.i 



-0.4 

3-1-0.3 

-0.2 



1.3:: 
2.4 ±0.5 
0.9 ±0.2 
1.2 ±0.3 

u-'-o.s 

-■-•'-0.5 

9.3 ±1.2 
3 1+^-1 



9.4 



-1-1.9 
1.6 



< 5 > (rg) 
(5) 



41^ 



fT2S 
-17 

7Atl1 
13.81^3 

9-1+29 
^-'--12 



NGC 4151. The mass determine d from re verberation map- 
ping is 4.57^0:47 X 10^ Mq ( .Bentz et al.|[2006 ). A stellar 
dynamical mass measurement obtamed using ground based 
I-band long-slit spectra and the Schwarzschild orbit superpo- 
sition code of Valluri et al. ( 2004 ) gave only an upper limit of 
^^ 4 X 10'^ Mq qOnicen et al.|| 2007p. However, a recent re- 



sult based on high resolution near infrared integral field spec- 
troscopic observations of the CO-band head with the Gemini 
North NIFS spectrograph, gives a revised stellar dynamical 
mass of (8.5 ± 3) x IO'^Mq (M. Valluri; private communi- 
cation) - nearly a factor of 2 greater than previous estimates. 
If this turns out to be right, Dcor for this source presented in 
Fig 2 should be closer to 5rg. 

As for the other two sources peaking at ^ 9rg, it is possi- 
ble that they were indeed caught at a time when their coro- 
nae were further than the median. It is interesting that RE 
J 1034+396 is the only persist ent AGN to ever sh ow a quasi- 
periodic oscillation (QPO; Gierlinski et al."2008). In the 
stellar-mass BH XTE J 1650-500, it has recently been shown 
that the onset of the more coherent QPOs are associated with 
the collapse of a corona from ^ lOrg during the transition 
from the ha rd-intermediate to the soft-intermediate state ( |Reis| 



etal. 



2013| lpl 



3. SIZE BASED ON MICROLENSING 

The possibility of using the "lens-like" effect of stars on dis- 
tant object is not new. However, its practical use in directly 
probing the inner regions around SMBHs was not established 
until the discovery of a flux-ratio anomaly in optical macroim- 
ages of the gravitationally lensed quasar QSO 2237+0305 
where the optical flux ratio between the bright est macroim- 
ages was statistically inconsistent with unity (Irwin et al. 
T989|j Although originally microlensing was not the only 
explanation for the flux anomaly problem in the optical (e.g. 
Dalai & Kochanek 2002 ), the systematic detection of an in- 
crease in the anomaly magnitude at higher energies strongly 
favoured microlens ing as the cause of these anomalies (e.g. 
Pooley et aLpOOT) !. As the effect of microlensing increases 
as the size scale of the lensed region decreases - as long as 
the source is smaller than the projection of the Einstein radius 
of the microlens into the plane (for a detailed review of mi- 
crolensing techniques see Wambsganss 2006 ) - this discovery 
confirmed the expected hierarchy in the characteristic sizes of 
the optical, UV and X-ray emission regions, with the latter 
often thought to be the most compact. 

Microlensing studies at optical/IR have o ften arrived at ac 



cretion disk sizes in the order of 10 r-g (e.g. Blackburne et al. 



The on ly other Q PO detected from a SMBH comes from a tidal disrup- 
1 source jReis et al. 20121 and is not persistently active. 



* Under a smooth gravitational potential, without the effect of micro- 
lensing, the brightest image-pairs are expected to be mirror-images of each 
other and have similar intensities iMetcalf & Zhao 2002 1. 
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Fig. 3. — Characteristic size of the X-ray emitting region versus mass (top) 
and characteristic luminosity (bottom). 



2006 1, and most remarkably, such studies have shown that 
this radius scales with BH mass in a manner consisten t with 
the expectations of thin-disk theory ( [Morgan et al.|2010 ). We 
are now at a stage where the size of the X-ray emitting corona 
in the innermost regi ons around these BHs are beginning to be 
measured (e.g..Morgan"et al.|2008l|Chartas et al.|2009[|Mos- 



m 






quera et al.|2013 and references therein). We list in lable 1 
(Bottom) the details of the 6 sources with a quantitatively de- 
rived half-light radius for the X-ray emitting region. These 
values were compiled fro m the various work l isted above and 
especially from Fig. 8 of Mosquera et al. ( 2013) . As these are 
all SMBHs, the accretion disk emission peaks in optical/UV, 
so the size of the X-ray emitting region therefore directly re- 
lates with the coronap] In compiling the various masses pre- 
sented in Table 2, we have searched the literature for the lat- 
est BH mass estimates based on the MgII(A2798A) mass- 
line width relation for all sources but HE 1104-1805, where 
CIV(A1549 A) was used inste ad. The values quoted are from 
the work of [Peng et al^ ( |2006l ) and |Sluseet aL] ( [20T2) i. 

In a similar manner to § 2, we have estimated a probability 
distribution for the size of the X-ray emitting region of each 
source in units of r^ based on 10^ randomly sampled values 
of mass and size (in cm) assuming either a Gaussian distri- 
bution with symmetric errors in logarithmic scale or, in the 
cases where such are not present, a uniform distribution over 
the range given in Table 1 . Figure 3 summarises the result for 
the characteristic sizes of the coronae. 

4. DISCUSSION 

The two techniques explored in this work effectively probe 
different aspects of the same emission region. From the rever- 
beration lags, we show in Fig 2 that 14/17 sources in our sam- 
pla^are at distances strongly clustered around 2-3 rg above 
the accretion disk. Taken at face value and assuming that the 
corona is radially symmetric about the central BH, this im- 
mediately suggest that a large fraction of Seyferts should be 
rapidly spinning (as a point of reference, a radius i?isco = 
3.15rg corresponds to a spin parameter a = 0.75cJ/ GAP), 
as is o ften found with detailed spectral analyses (Walto net al.| 
20T3). For a few sources, notably Mrk 841 and NGC 6860, 



' The use of half-light radius in microlensing studies is of particular impor- 
tance as these have been shown to be independent of the shape or the surface 
brightness profile of the source (e.g. Mortonson et al.|2005|for details). 

* 15/17 using the higher mass estimate of NGC 4151. 



the coronae appear to be at distances commensurate with radii 
for which there are no stable orbits, even for a maximally ro- 
tating BHj However, the effect of lag dilution c ould artifi- 
cially reduce the intrinsic lag (e.g. Wilkins & Fabian 2013} 
thus increasing the distance between the disk and coronEpl 
and possible errors in the BH mass could also change the ao- 
solute -Dcoi reported here. 

A number of other systematics uncertainties in the deriva- 
tion of soft-lags, including the unknown geometry of the 
corona or intrinsic flux variation (e.g.|Zoghbi et al.|2012[|Kara| 



et al.| 2013l could also increase the ra dius by a factor of a 
few. Indeed, for the case of NGC 4151, |Zoghbi et al.| ( |20T2] l 
showed that the lag-time differed by a factor of ^ 10 between 
high and low flux states. The result presented here conser- 
vatively shows the upper value (^ 2000 s), as found in the 
low-state of the source. If instead, we use the ^ 200 s lag, the 
distance between the disk and the corona would commensu- 
rate with a single gravitational radii. It is clear that for the ma- 
jority of the Seyferts reported here, the distance between the 
X-ray emitting region and the inner accretion disk is smaller 
than ^ 20r'g even if we allow for an order of magnitude sys- 
tematic shift in the derived values. The microlensing results 
further suggest that these coronae are highly compact, with 
upper limits on their characteristic sizes in the order of 20- 
30rg. 

The objects explored here span a large range in mass and 
Eddington fraction; however, in all cases our selection crite- 
ria (L > lO^^LEdd) was such that the systems are thought 
to be i n states where the accretion disk extends down to the 
ISCO ( jReis et al.||2013| and references therein);. Based on 
the presence of highly compact coronae, we postulate that the 
Seyferts and quasars presented here are in spectral states anal- 
ogous to the hard/intermediate and the soft/intermediate states 
of stellar mass BHBs. 

The compactness of the coronae in these sources strongly 
constrains models of Comptonisation, with magnetic recon- 
nection/flares (e.g. 'Merloni & Fabia n|200l] l or failed jet mod- 
els (e.g. Ghisellini et al. 2004) strongly favoured over mod- 



els which predict larger, extended coronae (e.g. Haardt & 
Maraschi_1991) . The small distance between the corona and 
the accretion disk further suggests that models where the en- 
ergy is dissipated in the innermost regions or possibly even 
within the plung ing regions, are fav oured. Such a scenario 
was presented by Hirose et al. (12004), where the corona is de- 
scribed as a region of smooth magnetic field lines increasing 
in strength toward small radii and being strongest for max- 
imally rotating BHs. More recently, global MHD simula- 
tions were used to show that the vast majority of the coro- 
nal emissio n comes fro m a relatively small volume of space 
(Schnittman et al.|2012[ ), consistent with observational results 
presented in this work. 

We also note that the base of the jet in M87 has recently 
been spatially resolved to be < 11 rg in size (Doeleman et al. 
|2012| l. If this is ass ociated with the coronae - as is the case m 
some models (e.g. |Falcke & Biermann|1999| l - then this too 
argues for a highly compact corona. 

A number of arguments have so far eluded discussion and 
shall be mentioned here only briefly. First, the spectral fits to 

' Again under the assumption that the corona is radially symmetric about 
the central BH . ^_^ 

'" Wilkins & Fabian|J2013} showed that this effect could increase Dcor 
in IH 0707-495 by up to 75%, and a factor of similar magnitude could also 
affect all other lags mentioned here. 



broad, relativistic lines profiles often require emissivities that 
are indicative of a very compact corona. This has been seen 
in both stellar-mass as well as SMBH and is often interpreted 
as the gravi tational bending of ligh t toward the inner regions 
of the disk ( |Miniutti & Fab ian 2004). We have chosen not to 
use this line of evidence as argument for a compact corona 
as it has recently been shown that these results could poten- 
tially be dependent on the prior assump tion that the corona is 
inherently compact ( Dauser et al.|2013] l. 

The presence of compact coronae are also often inferred by 
the existence of high frequency QPOs whose frequencies are 
similar to that expected from the Keplerian frequencies in the 
inner few gravitational radii from the black hole. QPOs are 
known to be stronger at higher energies which has prompted 
the association of the corona as the main instrument in mod- 
ulating the signals. Indeed, as mentioned above, QPOs have 
been seen in RE J1034H-396 and in the tidal disruption source 
Swift J164449.3H-573451 (Reis et al. 2012 ), and in both cases 
the QPOs have been associated with the coronae as the emis- 
sion from the disk peaks in the optical/UV. However, the 
drawback in inferring the properties of the corona from the 



QPO, is that to date there is no self-consistent model that 
can fully explain the observed behaviour in a wide variety of 
sources. 

In summary, current evidence - based on both imaging 
and timing techniques - strongly favours a compact corona 
located a few gravitational radii above the accretion disk. 
Whether the corona follows a slab, spherical, conical or 
patchy geometry, or whether it contains a purely thermal, non- 
thermal or hybrid population of electrons is still unknown. 
However, future emission models should ensure that the ma- 
jority of the energy dissipation occurs in a compact region 
close to the central black hole, if they are to be physical. Of 
course, both X-ray macrolensing and soft-lag techniques are 
still in their early days, and the application of these techniques 
in determining the size and location of coronal regions should 
be checked against a much larger sample in the near future. 
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